function W=Homosked_weight(data,theta)
%Calculates efficient weight matrix under the assumption of homoskedasticit
%(i.e. the 2SLS weight matrix in linear IV)
Z=data(:,3:end);
W=(Z'*Z)^-1;
end
